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I .  BACKGROUND 


AljjGa ^_jjAs/GaAs  MODFETs  have  become  increasingly  important  in  III-V 
semiconductor  device  circuit  development.^  Although  most  significant  at 
cryogenic  temperatures,  speed  and  power  advantages  over  GaAs  MESFET  techno- 
logics  are  evident  even  at  room  temperature.  In  applications  where 
epitaxial  growth  is  necessary,  such  as  GaAs  on  Si,  there  is  no  impediment 
to  the  use  of  these  devices.  As  GaAs  MODFETs  are  being  developed  for  a 
variety  of  high-speed  analog  ( >60  GH2)  and  digital  (200  MHz  clock  rate) 

•3  ij 

spaced-based  systems,-'’  an  understanding  of  the  radiation  degradation  of 
these  devices  is  of  particular  importance. 

In  an  attempt  to  understand  the  results  of  neutron  testing  of 
Aerospace  fabricated  MODFETs  we  have  developed  a  charge  control,  depletion 
layer  model  for  neutron-induced  threshold  voltage  shifts  in  AlGaAs/GaAs 
hetero junction  devices.  In  the  following  subsection  we  review  the  salient 
features  of  MODFET  band  structure.  The  subsequent  subsection  contains  a 
lengthy  summary  of  the  threshold  voltage  results  of  our  MODFET  strong  in¬ 
version  model,  which  has  been  discussed  elsewhere.^ 

Section  II  contains  the  details  of  the  model  used  to  describe  neutron 
degradation  in  these  structures,  and  how  the  neutron-induced  threshold 
voltage  shifts  are  derived.  In  Section  III,  experimental  results  for  the 
neutron- induced  threshold  voltage  shifts  are  compared  with  the  analytical 
results  derived  in  Section  II. 

A.  STRUCTURE 

The  band  structure  of  a  typical^  AiGaAs(n)/GaAs  heterojunction  with  a 
Schottky  gate,  under  bias,  is  shown  in  Fig.  1.  The  wide  bandgap  AlGaAs 
layer  of  thickness  (a+d)  lies  between  the  gate  and  the  narrow  bandgap  GaAs 
layer.  The  mismatch  between  the  bandgaps  causes  the  formation  of  a  two- 
dimensional  quantum  well  at  the  AlGaAs/GaAs  interface.  The  doped  AlGaAs 
layer  of  thickness  d  is  separated  from  the  two-dimensional  channel  by  a 
spacer  layer  of  thickness  a.  Doping  concentrations  in  the  AlGaAs(n)  layer 
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Fig.  1,  Pre-irradiation  Band  Diagram  of  a  Typical  MODFET  Structure 
with  Schottky  Contact,  4>m,  Under  Bias,  V  ,  AlGaAs  donors 
and  unintentional  acceptors  are  assumed  ®to  be  completely 
ionized.  Near  threshold  the  Fermi  level  (dashed  line)  lies 
below  the  bottom  of  the  2-D  well. 


1  fi  ^ 

are  on  tne  order  of  10  °  cm"-^  and  are  the  source  of  mobile  carriers. 

Because  the  2-D  channel  is  separated  from  the  highly  doped  AlGaAs(n)  region 
and  the  "unintentional"  acceptor  doping  concentrations  in  the  GaAs  are  low, 

10  1C  0 

typically  10  -  10 cm  for  molecular  beam  epitaxy  (MBE)  processed  de¬ 

vices,  enhanced  mobilities  are  achieved  (particularly  at  low  temperatures) 
due  primarily  to  reduced  impurity  scattering.^ 

B.  THRESHOLD  VOLTAGE  DEFINITION 


In  the  depletion  layer  approximation  under  the  conditions  depicted  in 
Fig.  1,  it  can  be  shown  by  application  of  Gauss's  law  across  the 
that  the  field  at  the  AlGaAs/GaAs  interface  is 


^i  =  (-^off  -  ^f 


(1) 
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where 


Voff  =  the  voltage  which  "annihilates"  the  channel  charge 

=  the  Fermi  level  relative  to  the  bottom  of  the  well 
V  =  applied  gate  voltage 

D 

d  =  the  doped  AlGaAs(n)  layer  thickness 
a  =  the  spacer  layer  thickness. 

The  term  is  given  by 

“off  =  -  ‘E,  -  v„  (2) 

where 

(ti^  =  the  Schottky  barrier  height 

AEg  =  the  conduction  band  offset 

V^2  =  voltage  drop  across  the  AlGaAs  layer. 

For  an  AlGaAs  donor  density  of  in  the  depletion  layer  approximation,  the 
voltage  01  up  across  tne  AlGaAs  layer  is 

Vai  =  -qNj2d^/2e  (minus  the  pinch-off  voltage,  Vp)  (3) 

where 

q  =  the  electronic  charge 

e  =  the  AlGaAs  permittivity  (assumed  to  be  identical  to  the  GaAs 
permittivity  in  what  follows). 

Substituting  Eqs .  (2)  and  (3)  into  Eq.  (1)  we  may  write 

Vg’  =  a)E2  f  E^^  (‘^) 
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where 


V  -  A 
g  ’rn 


In  the  neutron  degradation  analysis  that  follows,  the  Schottky  barrier 
height  and  the  conduction  band  offset  are  unaffected.  It  is  convenient, 
therefore,  to  consider  only  V 

O 

In  Appendix  A  we  show  that  the  Ferrtii  level,  E^,  depends  on  the  channel 
density,  n^,  in  the  following  way: 

=  Eq  +  (kT/q)ln(C2n3/(kT/q))  (5) 

where 

Eq  =  the  first  quantum  level  in  the  2-D  well 
kT/q  =  the  absolute  temperature  in  volts 

=  the  unintentional  acceptor  doping  density 
W  :  the  GaAs  depletion  layer  width. 


As  shown  in  Appendix  B  the  first  quantum  level  may  be  written  as 

Eo  =  Ci(n3.N3W)2''3  (6) 

In  MODFETs  we  take  threshold  to  occur  at  strong  inversion  when  the  channel 
carrier  density,  n„,  is  equal  to  the  local  acceptor  density,  N„,  times  the 
average  width  of  the  two-dimensional  channel,  Thus, 


•^th  ■  ^a^ave 


(7) 


r.iis  relationship  for  the  channel  density  at  threshold  is  consistent  with 
the  strong  inversion  definition  of  threshold  in  MOSFETs.  In  that  case, 
strong  inversion  occurs  when  the  density  of  electrons  in  the  channel 
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is  eai’ai  to  the  density  of  acceptors  in  the  bulk.  Our  definition  is  the 
two-dimensional  analog.  The  consequences  of  this  definition  are  discussed 
in  more  detail  elsewhere.^  In  Appendix  C  we  show  that  the  channel  density 
in  the  single  subband  triangular  well  approximation  at  threshold  is  related 
to  the  depletion  layer  width  in  the  following  way: 

"th  =  Co(N3)2/3/w1/3  (8) 


Substitution  of  Eqs.  (5),  (6),  (7),  and  (8)  into  Eq.  (4)  yields  V  '  at 

O 

threshold  in  terms  of  the  AlGaAs(n)  doping  density,  the  acceptor  doping 
density,  N^,  the  depletion  width,  W,  and  appropriate  geometrical  factors. 

As  in  the  case  for  MOSFETs,  the  depletion  width  is  related  to  the  band 
bending  in  the  bulk,  by 

bb 


[2(c/q)^BB/Na^ 


1/2 


(9) 


With  the  help  of  Fig.  1,  the  band  bending  is  shown  to  be: 


•^bulk 


-  E, 


(  10) 


I 


where 

Eg  ;  the  GaAs  bandgap 

((1,  ,  :  the  GaAs  bulk  potential  [-  (kT/q)ln(N- /n^  )  ]  (n^  is  the  GaAs 

bulk  A  _'3  ®  ^  ^ 

intrinsic  carrier  density  -  10  cm  at  room  temperature) 

This  result  for  the  band  bending  in  MODFETs  is  slightly  different  than  that 
for  MOSFETs  and  is  discussed  elsewhere.^  With  the  use  of  Eqs.  (9)  and 
(10),  Vg'  at  threshold,  which  we  will  write  as  depends  only  on  the 

doping  densities  and  the  geometry. 
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Application  of  Gauss's  law  from  the  interface  to  the  edge  of  the 
depletion  layer  yields  the  field  at  the  interface  in  terms  of  GaAs  doping 
only.  Inspection  of  Eq.  (4)  shows  that  the  contribution  to  the  threshold 
voltage  by  the  AlGaAs  material  parameters  and  the  GaAs  material  parameters 
may  be  separated  and  examined  independently.  This  property  of  the  thres¬ 
hold  voltage  will  be  exploited  when  the  neutron-induced  threshold  voltage 
results  are  discussed. 


10 


II.  NEUTRON-INDUCED  THRESHOLD  VOLTAGE  SHIFTS 


After  neutron  irradiation  and  at  threshold,  Eq.  (4)  may  be  rewritten 
to  reflect  the  post- irradiation  value  of  Vth’  thus 

''th"  =  '^Al"  *  *  a)Ei"  Ef"  (11) 

where  double  primes  (")  indicate  the  post-irradiation  values.  The  differ¬ 
ence  between  the  pre-irradiation  value  ,  and  is  the  neutron- 

induced  threshold  voltage  shift,  or 

^''th  =  ^^Al  (d  a)AE^  +  AEf.  (12) 

In  the  following  subsection  the  model  which  describes  the  redistribu¬ 
tion  of  charges  due  to  neutron  bombardment  is  described.  The  subsequent 
subsection  contains  a  discussion  of  threshold  voltage  shifts  due  to  each  of 
the  terms  in  Eq.  ( 12) . 

A.  NEUTRON- INDUCED  TRAP  LEVELS  AND  CHARGE  DENSITIES 

It  is  well  known  that  high  energy  neutron  bombardment  of  semiconduc- 

8  Q 

tors  introduces  deep  traps  in  the  bandgap  of  the  material.  In  lieu  of 
any  detailed  knowledge  of  the  exact  location  of  these  traps,  we  assume  that 
the  neutron-induced  trap  may  be  described  as  an  effective  amphoteric  trap 
located  at  midgap.  This  assumption  reflects  exper imeiital  results  which 
show  that  neutron  bombardment  of  both  p  and  n  type  material  results  in  the 
material  becoming  more  intrinsic. Making  this  assumption  for  both  the 
AlGaAs  and  the  GaAs,  the  MODFET  band  diagram  after  neutron  irradiation  may 
be  depicted  as  shown  in  Fig.  2.  Amphoteric  traps  located  below  the  Fermi 
level  charge  negatively  whereas  traps  above  the  Fermi  level  charge  posi¬ 
tively.  As  a  result,  the  voltage  drop  across  the  AlGaAs,  the  field  at 

the  interface,  E^,  and  the  Fermi  level,  E^-,  at  threshold  will  change  with 
neutron  irradiation. 
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Fig.  2,  Post-irradiation  Band  Diagram  of  a  Typical  MODFET 

Structure.  Effective  amphoteric  midgap  traps  charge 
negatively  when  below  the  Fermi  level  and  charge 
positively  when  above  the  Fermi  level.  At  threshold, 
some  midgap  traps  in  the  GaAs  lie  below  the  Fermi 
level  and  some  lie  above.  W_  is  the  extent  of  the 
region  in  the  GaAs  in  which  the  traps  are  below  the 
Fermi  level.  is  the  extent  of  the  region  in  which 
the  traps  lie  above  the  Fermi  level. 


We  assume  that  the  decrease  in  the  net  charge  in  the  AlGaAs  layer  due 
to  charging  of  the  neutron-induced  traps  may  be  effectively  characterized 
as  the  removal  of  carriers.  To  first  order  carrier  removal  has  been  shown 
to  be  linear  with  neutron  fluence.^^  Therefore,  the  net  charge  in  the 
doped  AlGaAs  layer  after  neutron  irradiation  may  be  written  as 


Nd'  =  Nd  -  ad* 


(13) 
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where 


:  the  carrier  removal  rate  in  the  AlGaAs 
4)  :  the  neutron  fluence. 


In  a  more  detailed  treatment  of  the  charging  of  eunphoteric  traps,  the 
carrier  removal  rate,  a^j,  is  calculated  by  accounting  for  the  population 
statistics  of  these  traps  and  averaging  the  charge  over  the  extent  of  the 
AlGaAs.  As  we  intend  to  treat  the  carrier  removal  rate  as  an  experiment¬ 
ally  determined  parameter,  we  will  dispense  with  such  detail . 

At  threshold  it  is  possible  that  some  of  the  amphoteric  traps  in  the 
GaAs  layer  lie  below  the  Fermi  level  and  that  some  lie  above  the  Fermi 
level.  We  assume  that  in  each  region  the  net  charge  is  linearly  dependent 
on  the  neutron  fluence,  to  first  order.  In  keeping  with  the  depletion 
layer  approximation  the  net  charge  density  in  the  GaAs  is  then 


N-  •••  a  <ti*,  (d+a)<z<(W  +d+a) 

d  - 

(1A) 

-  3^4;  (W_+d+a)<z< (W^+W_+d+a) 

(15) 

where 

a_  =  is  the  effective  acceptor  introduction  rate  in  the  GaAs 

a_^  =  is  the  effective  carrier  removal  rate  in  the  GaAs 

and 

W_  =  the  region  in  the  GaAs  in  which  the  traps  lie  below  the  Fermi 
level  and  therefore  charge  negatively 

=  the  region  in  which  the  traps  lie  above  the  Fermi  level  and 
therefore  charge  positively. 


13 


B.  TRESHOLD  VOLTAGE  SHIFT 

Given  the  post-irradiation  charge  density  in  the  AlGaAs  layer,  N^j', 
the  voltage  drop  across  this  layer  is 

Vai’  =  -qN^’d^/Be  (16) 

by  analogy  with  Eq.  (3).  The  difference  between  Eqs.  (16)  and  (3)  yields 
the  threshold  shift  due  to  the  AlGaAs  layer: 

AVai  =  qa^4>d2/2  .  (17) 

It  should  be  pointed  out  that  we  have  ignored  charging  of  the  traps  in 
the  spacer  layer  to  derive  the  result  of  Eq.  (17).  The  contribution  to  the 
threshold  voltage  shift  due  to  the  spacer  layer  is  about  7J  of  that  deter¬ 
mined  from  Eq.  (17).  It  will  be  shown  later  that  the  AlGaAs  contribution 
to  the  threshold  voltage  shift  is  less  than  of  the  total,  therefore  the 
spacer  layer  contribution  is  negligible  (<<1yt).  The  possibility  that  some 
traps  in  the  AlGaAs  might  lie  above  midgap,  and  charge  positively,  has  also 
been  ignored.  This  phenomenon  would  tend  to  lessen  the  contribution  to  the 
threshold  voltage  shift  from  the  AlGaAs  layer  even  more  because  the  net 
neutron- induced  charge  across  the  AlGaAs  layer  would  be  smaller.  Because 
the  carrier  removal  rate,  a^,  is  determined  by  electrical  measurements, 
these  factors  (spacer  layer  charging  and  positive  neutron- induced  charges 
in  the  AlGaAs  layer)  are  implicitly  accounted  for. 

Application  of  Gauss's  law  from  the  interface  into  the  bulk  of  the 
GaAs  with  the  post-irradiation  charge  densities  yields 

E^"  =  (q/e)[(N3  +  a_4>)W.  +  (Ng  -  a^^)W^]  (18) 

for  the  field  at  the  interface.  As  written  in  Eq.  (18),  the  contribution 
to  the  field  at  the  interface  due  to  the  channel  charge  has  been  ignored. 

At  threshold  the  contribution  to  the  field  by  the  channel  charge  is  of  the 
order  z^ve^^^-'  channel  charge  width,  is  on  the  order  of  10  nm 
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and  the  depletion  layer  width,  W_,  is  on  the  order  of  1  urn.  Therefore, 
contributions  of  order  dropped.  The  pre- irradiation 

interface  field  may  be  calculated  to  give 

Ej  =  (q/ONgW.  (19) 

The  difference  between  Eqs.  (19)  and  (18)  times  the  AlGaAs  layer  thickness, 
(d+a),  yields  the  neutron-induced  threshold  voltage  shift  due  to  the  redis¬ 
tribution  of  charge  in  the  GaAs  layer. 

Deep  in  the  GaAs  layer,  after  neutron  irradiation,  the  charge  density 
is  given  by  the  second  term  in  Eq.  ( 18)  as 

0^  =  (20) 

At  high  fluences  the  neutron-induced  positive  charge  buildup  reduces  the 
net  charge  in  the  layer.  Therefore,  the  field  contribution  of  the 
charges  in  the  layer  relative  to  the  the  contribution  due  to  the  bulJdup 
of  negative  charges  in  the  W_  is  reduced.  Later,  we  show  that  the 
threshold  voltage  shift  in  MODFETs  is  not  appreciable  for  low  fluences 
(<10  ^  cm*  ) .  Therefore,  we  shall  ignore  the  contribution  to  the  field  at 
the  interface  due  to  the  charges  in  the  region.  The  change  in  the  field 
at  the  interface  may  be  written  as 

AE^  =  (q/e)[(Ng  a_*)W_  -  N^WJ.  (21) 

Because  the  effective  acceptor  introduction  rate,  a_,  will  be  treated  as  an 
experimentally  determined  parameter  as  long  as  the  net  charge  buildup  in 
the  GaAs  is  negative,  Eq.  (21)  may  be  used  to  describe  the  change  in  the 
field  at  the  interface  due  to  the  neutron  induced  charge  redistribution. 

After  irradiation  the  Fermi  level,  E^.",  may  be  determined  from  Eq.  (5) 
by  evaluating  the  post-irradiation  threshold  channel  charge  density,  n^j^', 
using  the  post-irradiation  effective  charge  density,  N_,  and  substituting 
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into  Eq,  (5).  Subtracting  the  pre-irradiation  Fermi  level  from  the  post- 
irradiation  level  and  ignoring  terms  of  order  2j^yg/W_  yields 

AEf>  =  (kT/3q)ln[(1  +  a.o/Ng^)2w/W_]  (22) 

for  the  Fermi  level  contribution  to  the  neutron- induced  threshold  voltage 
shift. 

In  order  to  fully  descrioe  the  neutron- induced  threshold  voltage 
shift,  it  remains  to  determine  W_  in  terms  of  GaAs  parameters.  Application 
of  Gauss's  law  from  the  interface  into  the  bulk  of  the  GaAs  yields  the 
following  for  the  change  in  potential  from  the  interface  to  the  edge  of  the 
W_  region: 


V^-V_  =  (q/£)({N3-.a_4>)W_2/2  +  {N^-a^^)W^W_  ] .  (23) 

Ignoring  the  net  charge  in  the  region  compared  to  the  charge  in  the  W_ 
region  as  discussed  earlier,  the  width  of  the  W_  region  may  be  determined 
as 


W.  =  [2c(V^-V_)/q(N3  +  a.«)]’/2^  (24) 

As  seen  in  Fig.  2,  the  potential  change  across  the  W_  region  may  be  written  as 

Vi-V_  =  Eg/2  -  Ef".  (25) 

Because  the  neutron- induced  change  in  the  Fermi  level  is  small  compared  to 
E  /2,  this  may  be  approximated  by 

D 

Vi-V.  =  Eg/2  -  Ef.  (26) 

Substitution  of  Eqs.  (2M)  and  (26)  into  Eq.  (21)  and  using  the  earlier 
results  for  the  depletion  layer,  W,  yields  the  neutron-induced  threshold 
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voltage  shift  as  a  function  of  pre-irradiation  material  parameters,  geome¬ 
try,  the  effective  acceptor  introduction  rate,  the  carrier  removal  rate, 
and  the  neutron  fluence. 

C.  LOW  FLUENCE  EXPANSION 

Due  to  the  cumoersome  nonlinear  dependence  of  the  threshold  voltage 
shift  on  neutron  fluence,  a  discussion  of  the  relative  importance  of  each 
term  in  Eq.  (21)  in  the  low  fluence  limit  is  warranted.  Expanding  the 
neutron-induced  threshold  voltage  in  a  Maclaurin  series  for  small  fluence, 
such  that  a_4i/N^  <v  1,  and  keeping  first  order  terms  only,  the  threshold 
voltage  shift  may  be  written  as 

'“th  =  127) 

where  is  given  by  Eq.  (3)  and  depends  on  GaAs  material  parameters 

in  the  following  way: 

^^"Ca  ■  I(d  +  a)  dAE|/d4>f  dAE^/di^]  iji  =  0  (28) 


where 


[diE^/d«]«  :  0  =  (q/e)a_[W_]4)  =  0  (29) 

[dAE|./dt]ct  =  0  =  (5kT/6q){a_/N^).  (30) 

For  typical  MODFET  geometries  (<  500  A  AlGaAs),  assuming  that  the 
effective  acceptor  introduction  rate  and  the  carrier  removal  rate  are  of 
the  same  order  (to  be  justified  in  Section  III),  and  that  the  unintentional 
acceptor  doping  is  much  less  than  the  AlGaAs(n)  doping  (a  ratio  of  -10~^  is 
typical),  the  AlGaAs  contribution  to  the  neutron-induced  threshold  voltage 
shift  is  of  the  order  of  This  low  fluence  expansion  applies  to 

fluences  of  less  than  5  *  10^^  cm"^.  However,  the  trend  of  GaAs  dominance 
of  the  neutron  degradation  continues  to  fluences  well  above  1  »  10^^  cm"^ 
as  shown  in  Section  III. 
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III.  EXPERIMENTAL  COMPARISON 


MODFETs  were  fabricated  on  commercially-obtained,  MBE  grown,  epitaxial 
material.  The  epilayers  consisted  of  1  um  unintentionally  doped  GaAs, 

20  A  unintentionally  doped  AIq  25^^0  75^^,  500  A  n-AlQ  25^®0  75^®  " 

1  X  10^®  cm“^  Si),  and  200  A  n-GaAs  (n  =  1  x  10^®  cm“^  Si).  The  aluminum 
mole  fraction  of  x  =  0.25  was  chosen  to  minimize  complications  from  high 
concentrations  of  pre-irradiation  deep  levels  characteristic  of  higher  mole 
fraction  alloys.  A  standard  mesa  isolation  gate  process  was  used  to  define 
the  transistors.  Figure  3  shows  a  stylized  cross  section  of  the  devices 
fabricated  for  this  study. 

Ili 

Using  an  unintentional  acceptor  doping  density,  N^,  of  5  *  10  cm“^, 

AlGaAs(n)  doping,  N^^,  of  1  x  10^®  cm~^,  a  post-processed  doped  AlGaAs  layer 
thickness,  d,  of  i400  A,  and  a  spacer  layer  thickness,  a,  of  20  A  our  thres¬ 
hold  voltage  analysis,  yields  a  value  of  -.620  volts  for  the  pre- 
irradiation  threshold  voltage.  This  is  in  reasonable  agreement  with  our 
measured  values  which  range  from  -.41  to  -.73  volts.  Variation  in  the 
measured  values  is  a  result  of  variations  in  the  thickness  of  the  AlGaAs 
arising  from  non-uniform  processing.  The  value  used  for  analysis, 

420  A,  is  a  nominal  value  and  gives  rise  to  the  discrepancy  between  our 
calculated  value  and  the  measured  values. 

The  devices  were  exposed  to  neutron  fluences  ranging  from  5  *  10^^ 

-2  m  -2 

cm  to  1  >  10  cm  in  the  SPR-III  nuclear  reactor  at  Sandia  National 
Laboratories.  Shown  in  Table  1  are  the  measured  threshold  voltage  shifts 
for  the  Aerospace  fabricated  MODFETs.  The  nominal  device  parameters  are 
those  cited  above.  In  Fig.  4  the  data  for  the  measured  threshold  voltage 
shift  from  Table  1  is  plotted  versus  neutron  fluence.  Also  shown  in  Fig.  4 
are  the  modeling  results  for  the  threshold  voltage  change,  assuming  an 
AlGaAs  carrier  removal  rate,  a^,  of  10  cm"^  and  a  GaAs  acceptor  introduc¬ 
tion  rate,  a_,  of  3  cm"\  The  solid  curve  in  Fig.  4  is  the  modeling  result 
of  the  full  model,  channel  charge  and  AlGaAs  layer  contribution  included, 
while  the  dashed  curve  is  the  result  obtained  by  making  the  simplifying 
assumptions  discussed  earlier  and  ignoring  the  AlGaAs  contribution. 
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Fig.  3.  Styli2ed  Version  of  the  MODFET  Structure  Used  in 
This  Study.  Low  mole  fraction  AlGaAs  was  used 
to  avoid  deep  pre-irradiation  electron  traps. 


Table  1.  Measured  Threshold  Voltage  Shifts 
of  Aerospace  Fabricated  MODFETS 


Neutron 

Threshold 

Fluence 

Shifts 

( Volts) 

9.2E  + 

0. 1864 

±  0.0214 

i|.9E  + 

0.1100 

±  0.0100 

2.2E  +  14 

0.0540 

±  0.0150 

1.0E  +  14 

0.0319 

±  0.0155 

5. IE  +  13 

0.0184 

±  0.0059 
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NEUTRON  FLUENCE  (cm 


Fig.  4.  Change  in  Threshold  Voltage  vs  Neutron 
Fluence  for  the  Model  Parameters  Shown 

For  these  values  of  a^j  and  a_,  the  agreement  between  the  model  and  the 

experiment  is  very  good  over  the  whole  range.  The  value  of  was  chosen 

from  the  literature^^  and  verified  by  Hall  measurement,  then  a_  was  used  to 

1  ii  -? 

fit  the  data.  Only  at  the  highest  fluence  levels,  -  5  »  10  cm  ,  do  the 
simplified  model  results  diverge  significantly  fron  the  full  model  results. 
The  modeling  results  for  the  threshold  voltage  using  these  parameters  are 
shown  in  Fig.  5.  The  increase  in  the  threshold  voltage  with  neutron 
fluence  is  very  similar  to  the  increase  in  threshold  voltage  due  to 

c 

increased  acceptor  doping,  discussed  elsewhere.-^ 

Shown  in  Fig.  6.  are  the  contributions  of  each  term  described  in 
Eq.  (11).  The  threshold  voltage  shift  data  and  modeling  results  of  Fig.  4 
are  included.  The  change  in  the  Fermi  level,  AE^,  and  the  change  in  the 
field  at  the  interface  to  the  redistribution  of  charge  in  the  GaAs,  E^ , 
dominate  and  are  almost  equal,  whereas  the  AlGaAs  contribution  is  about 
2.5/t  over  the  whole  range  of  neutron  fluence. 
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Fig.  5.  Threshold  Voltage  vs  Neutron  Fluence  Using  the  Model 
Parameters  Used  to  Fit  the  Threshold  Shift  Data 


MODEL  PARAMETERS 

CRR.  AIR  10.'3.  DOPING  D/A:  lE  +  18  5E*14 
AlGaAs 'SPACER  LAYER  THICKNESS  420:20  A 


Fig.  6.  Contributions  to  the  Threshold  Voltage 
Shift  vs  Neutron  Fluence 
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IV.  DISCUSSION 


In  Fig.  7  we  have  shown  a  diagram  of  the  conduction  band  of  the 
AlGaAs/  CaAs  interface  of  a  typical  MODFET  structure  at  threshold  as  the 
neutron  fluence  varies.  Also  shown  are  the  positions  of  the  Fermi  energy 
and  the  first  energy  level  in  the  2-D  well  as  the  fluence  increases  from 
some  initial  value,  (i),  to  some  final  value,  (f).  As  the  neutron  fluence 
increases,  the  Fermi  level  shifts  upward  closer  to  the  bottom  of  the  well 
and  the  first  energy  level  moves  up  away  from  the  bottom  of  the  well.  The 
difference  between  the  Fermi  level  and  the  first  quantum  level  determines 
the  amount  of  charge  in  the  channel.  In  keeping  with  the  definition  of 
threshold,  the  channel  charge  must  increase  as  the  effective  charge  density 
near  the  interface  increases  with  neutron  fluence  as  the  midgap  traps  near 
the  interface  charge  negatively. 

As  the  charge  in  the  channel  increases,  the  field  at  the  interface 
increases  and  the  2-D  well  width  decreases,  consistent  with  the  increase  in 
the  energy  level.  At  the  same  time  the  depletion  width  is  decreasing. 

Since  the  depletion  width  contribution  to  the  threshold  voltage  shift 
depends  on  the  product  of  the  depletion  width  and  the  effective  acceptor 
density,  an  increasing  depletion  contribution  means  that  the  increase  in 
the  effective  acceptor  density  must  dominate. 

For  the  device  parameters  cited  above,  we  have  plotted  the  Fermi 
energy  and  first  quantum  well  energy  level,  the  change  in  the  Fermi  energy 
and  the  change  in  the  first  quantum  level,  and  the  depletion  width  as  a 
function  of  neutron  fluence  in  Figs.  8,  9,  and  10  respectively.  Figure  8 
shows  that  both  the  Fermi  level  and  the  first  quantum  level  increase  with 
neutron  fluence.  The  Fermi  level  increases  more  than  the  first  energy 
level,  as  shown  in  Fig.  9,  increasing  the  amount  of  charge  in  the  chan¬ 
nel.  Shown  in  Fig.  10  is  a  three-fold  decrease  in  the  depletion  width  from 
a  fluence  of  1  x  10^^  to  1  x  10^^  cm*^. 
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RFLATIVE  TO  Ttl(  ROFinM 


nFPlFlinN  WIDTH  (cm)  CHANGE  IN  ENERGY  LEVEL  (Volts) 


UU  g;=r=— .  .  ■■■■  L_! _ i _ I _ ^ _ L_J 

1E^13  1E  +  14  1E  +  15 

NEUTRON  FLUENCE  (cm"^) 


9.  Change  in  the  Quantum  Well  Level  and  the  Fermi 
Level  vs  Neutron  Fluence 


NEUTRON  FLUENCE  (cm'^) 


Fig.  10.  Depletion  Layer  Width  vs  Neutron  Fluence 
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V.  SUMMARY 


We  have  developed  and  presented  a  finite  temperature,  charge  control 
model  for  the  threshold  voltage  of  modulation-doped  field  effect  transis¬ 
tors  which  includes  degradation  due  to  neutrons.  The  model  has  been  ap¬ 
plied  to  Aerospace  fabricated  AlGaAs/GaAs  MODFETs  which  have  been  exposed 
to  neutron  fluences  from  5  *  10^^  to  1  >  10^^  cm~^.  Our  analysis  shows 
that  the  degradation  of  the  threshold  voltage  is  due  primarily  to  a  redis¬ 
tribution  of  charge  in  the  GaAs  rather  than  degradation  in  the  AlGaAs. 

This  implies  that  the  GaAs  doping  density  in  the  depletion  layer  is  an 
important  factor  in  determining  neutron  hardness  of  MODFETs.  Assuming  a 
constant  acceptor  introduction  rate,  higher  doping  in  the  GaAs  will  make 
these  structures  harder  to  neutrons. 
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APPENDIX  A  -  (Fermi  Level  and  n^) 


In  a  2-D  electron  gas,  the  energy  of  an  electron  in  the  ith  band  is 
given  by: 

Ei  =  Ejq  +  (A1) 

where 

M  =  the  Planck  constant  divided  by  2n 

k  =  the  longitudinal  wavevector 

=  the  longitudinal  effective  mass  of  the  carriers 
E^q  =  the  k  =  0  energy  level  of  the  quantum  well. 


The  concentration  of  carriers  is  calculated  by  evaluating  the  following: 


ng  =  ggd/A)  r  J  g{Epf(Ei)dEi 


i  E. 
10 


where 


gg  =  the  spin  degeneracy  of  the  carriers 
g{E^)  =  density  of  states 

A  =  device  area 

f(E^)  =  Fermi-Dirac  distribution  function 

I  =  the  sum  over  2-D  energy  levels, 

i 

Equation  (A2)  may  be  evaluated  by  considering  that: 

g{Ej)  =  g(k)dk  (A3a) 
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where 


g(k)dk  =  dAjj/Ajj 


(A3b) 


and 


=  Tik^:the  area  in  reciprocal  space  [=(2'h)^/A].  (A3c) 

Using  Eqs.  (A1)  and  (A3),  Eq.  (A2)  takes  the  following  form: 

a> 

"s  =  (gsmi/2ith^)  r  J  {1  +  exp[(E  -  Ef)/kT])  dE.  (Ai4) 

i  E. 

10 

Effecting  the  integration  in  Eq.  (A4),  assuming  only  one  energy  level,  E^, 
and  using  the  spin  degeneracy  factor,  gg,  for  electrons  (gg  :  2)  the  Fermi 
energy  may  be  calculated: 

Ef.  =  Eq  +  (kT)ln[exp(iilfl2n^/(n,^kT))  -  1)].  (A5) 

Near-threshold  ng  is  small  and  the  exponential  may  be  expanded.  Conversion 
to  potential  units  then  results  in  Eq.  (5)  in  the  text. 
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APPENDIX  B  -  (Eq  and  n^) 

The  first  quantum  level  may  be  calculated  in  the  triangular  well 
approximation  using  a  normalized  variational  envelope  wave  function: 

4i  =  (Bl) 

o 

where  b  is  a  variational  parameter  to  be  determined.  The  Hamiltonian  for 
the  triangular  well  is; 

H  =  -(h^/2m^)d^/dz^  +  qV(z)  (B2a) 

where 


V{z) 

:  E^z:  2  >  0 

(B2b) 

V(z)  = 

infinity:  z  <  0 

(B2c) 

h  =  the  Planck  constant  divided  by  2  n 

m^  =  longitudinal  effective  mass  of  the  carriers. 

Minimi2ing  the  energy,  with  respect  to  b  yields: 

b  =  [aEm^q/h^J^^^  (B3) 

where 

E  =  the  field  in  the  potential  well. 


The  first  energy  level  may  be  calculated  by  evaluating  the  following: 


Eo  =  «>ol“l^o> 

which  yields: 

Eq  =  (9H2/8m^q)(i<Emjq/h2)2/3.  (BUb) 

Application  of  Gauss's  law  from  the  interface  into  the  bulk  of  the  struc¬ 
ture  shown  in  Fig.  1  gives  the  field  at  the  interface: 

E  =  (q/eXng  +  N^W).  (B5) 

Substitution  of  Eq.  (B5)  into  Eq.  (B^jb)  yields  Eq.  (6)  in  the  text  where. 

=  (9M^/8mj^q)(4qm^q^/h^e 


APPENDIX  C  -  N^,  and  W) 


Assuming  the  criteria  for  threshold  as  stated  in  the  text: 

”th  "  **^a^ave 

where 


(Cl) 


N„  =  the  acceptor  density 

a 

^ave  ~  average  position  of  carriers  in  the  2-D  well. 

The  average  position  may  be  calculated  in  the  triangular  well  approx imat ion 
using  the  variational  wave  function  from  Appendix  B.  Thus: 

"ave  =  (C2a) 

which  reduces  to: 

=  3/6  (C2b) 

where  the  variational  parameter,  b,  is  given  in  Appendix  B.  Combining 
Eqs.  (Cl)  to  (C3)  and  substituting  for  b  yields: 

Zave  =  (C3) 

where 

E  =  is  the  field  in  the  potential  well. 
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At  threshold  the  field  is  dominated  by  the  depletion  layer  charge  and  may 
be  approximated  as: 


E  =  qNjjW/e 


(Cii) 


where 


W  =  the  depletion  layer  width. 


Substitution  of  Eqs.  (C4)  and  {C3)  into  Eq.  (Cl)  yields  Eq.  (8)  used  in  the 
text . 
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LABORATORY  OPE RAT IONS 


The  Aerospace  Corporation  functions  as  an  ''architect'englneer  '  for 
national  security  projects,  specializing  in  advranced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  Investigations  is  the  technical  staff's  wide-^ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-speclf Ic  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -fleld-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory;  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communlcactons; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  (oromunlcatlon  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory;  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  Che  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


